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ABSTRACT 



We present Extreme Ultraviolet Explorer spectroscopy and photometry of the 



nearby F8 V star HD 35850 (HR 1817). The EUVE short- wavelength 75-175 A and 



medium-wavelength 160-365 A spectra reveal 28 emission lines from Fe IX and Fe XV 
to Fe XXIV. The Fe XXI AA102, 129 ratio yields an upper limit for the coronal electron 
density, logn e < 11.6 cm . The EUVE SW spectrum shows a small but clearly 
detectable continuum. The 75-150 A line-to-continuum ratio indicates approximately 
solar Fe abundances, with 0.8 < Z < 1.6 (90% confidence interval). Upper limits have 
been derived for a dozen high-emissivity Fe X through Fe XIV lines. The resulting EM 

ST 

7.4. Over the course of the 1-week observation, large-amplitude, long-duration flares 
were not seen in the EUVE Deep Survey light curve, though the light curve does show 
signs of persistent, low-level flaring, and possible rotational modulation. 



The EUVE spectra have been compared with non-simultaneous ASCA SIS spectra 
of HD 35850 obtained in 1995. The SPEX DEM analysis of the SIS spectrum indicate 
the same temperature distribution as the EUVE DEM analysis. However, the SIS 
spectra suggest sub-solar abundances, 0.34 < Z < 0.81. Although some of the 
discrepancy may be the result of incomplete X-ray line lists, we cannot explain the 
disagreement between the EUVE line-to-continuum ratio and the ^.SC/l-derived Fe 
abundance. 

The X-ray surface flux on HD 35850 is comparable to that of cooler dwarfs of 
comparable age and rotation like EK Draconis (GO V) and AB Doradus (Kl V). Given 
its youth (i ~ 100 Myr), its rapid rotation (vs'mi ~ 50 km s -1 ), and its high X-ray 
activity (Lx ~ 1.5 x 10 30 ergs s _1 ), HD 35850 may represent an activity extremum 
for single, main-sequence F-type stars. The variability and EM distribution can be 
reconstructed using the continuous flaring model of Giidel provided that the flare 
distribution has a power-law index a ~ 1.8. Similar results obtained for other young 
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solar analogs suggest that continuous flaring is a viable coronal heating mechanism on 
rapidly rotating, late-type, main-sequence stars. 

Subject headings: stars: individual (HD 35850) - stars: coronae - stars: late-type - 
stars: rotation - X-Rays: stars - ultraviolet: stars 

1. Introduction 

By analogy to the Sun, magnetic heating in the chromospheres, transition regions, and 
coronae of late-type stars is observed as cooling in the Balmer lines, the Ca II lines, and in 
UV, FUV, EUV, and soft X-ray continuum and line emission. Results from the Mt. Wilson 
Ca II H and K survey of single lower main-sequence field stars has provided compelling evidence 
of Solar-like chromospheric active regions, activity cycles, differential rotation, and occasional 
Maunder minima in mid-F to early-M stars (e.g., Noyes et al. 1984; Baliunas et al. 1996; Gray & 
Baliunas 1997; Baliunas et al. 1995). 

Coronal heating and magnetic dynamos on the Sun and in stars have been discussed 
extensively in the literature (cf. Haisch &: Schmitt 1996). The various flares, microflares, 
nanoflares, and brightenings observed on the Sun are thought to be manifestations of current 
sheet reconnection and much of the coronal heating debate is focussed on whether or not these 
transient events can account for the radiative output of the solar corona (e.g., Oreshina & Somov 
1998). While both continuous (e.g., MHD waves) and stochastic (nanoflaring) processes probably 
contribute to magnetic coronal heating, the dominant process has yet to be identified. 

Empirically, coronal emission on active, cool, main-sequence stars is similar to flare emission 
from solar active regions. For example, Benz & Giidel (1994) find a similar correlation between 
X-ray luminosity and non-thermal continuum radio luminosity for solar flares and active stars. 
On the Sun, mildly relativistic electrons accelerated along coronal magnetic field lines produce 
gyrosynchrotron radio emission. The presence of persistent but variable nonthermal radio emission 
from cool stars indicates a continual replenishment of the relativistic electron population while a 
correlation between coronal and radio emissions suggests a causal relationship between electron 
acceleration in magnetically confined loops and coronal heating. In the 50,000-200,000 K regime, 
further evidence of continual flaring comes from the broad emission-line components seen in 
high-resolution UV spectra of active dwarfs and RS CVn binaries. Because the broad-component 
line widths are 2-4 times the thermal width and because the broad components are blue-shifted 
relative to the narrow components, Wood et al. (1996) suggest that the broad components are 
produced by transition-region explosive events. 

Two well-studied examples of extreme main-sequence magnetic activity are the Pleiades-age, 
rapid rotators AB Doradus (K0 V) and EK Draconis (GO V). Spectro-polarimetric monitoring 
of AB Dor shows long-lived, cool magnetic spots and latitudinal differential rotation (Donati & 
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Collier Cameron 1997). ROSAT PSPC soft X-ray photometry of AB Dor (Kiirster et al. 1997) 
shows rotationally modulated flare and quiescent emission with the same phase and period (12.4 h) 
as the photospheric spots. ASCA and EUVE spectra of AB Dor and EK Dra indicate coronal 
emission- measure distributions, Em(T), with peaks at 5-8 MK and 20-30 MK (Mewe et al. 1996; 
Giidel et al. 1997). 

The 5-8 MK component has been observed in other moderately active coronal sources, most 
notably Capella (Brickhouse, Raymond, & Smith 1995). As has been pointed out by Gehrels & 
Williams (1993), optically thin plasma tends to accumulate at temperatures where the cooling 
curve fi(T) has a positive slope. This is a manifestation of the Parker (1953) Sz Field (1965) 
instability: when > 0, plasma cools more efficiently at slightly higher temperatures and 

less efficiently at slightly lower temperatures. At temperatures where < 0, plasma quickly 

cools to lower temperatures. For solar- abundance plasmas, rf ^P > in the ranges 0.7-1.2 MK, 
4.5-7 MK, and above 21 MK. Thus, the shape of the cooling curve explains the dip from 2-4 MK 
and from 9-20 MK seen in many EM distributions and provides some physical justification for 
two- and three-temperature coronal models. The temperature of the hot component and the 
relative amount of hot emission measure must reflect a balance between magnetic heating and 
radiative cooling. For stars with significant emission measure above 20 MK, some mechanism 
must be super-heating the coronal plasma at fairly regular intervals. 

For three active, solar-mass stars observed with the ASCA SIS (EK Dra, HN Pegasii, and 
k 1 Ceti in order of decreasing activity), Giidel (1997) has modeled Em(T) by assuming X-rays are 
produced by an ensemble of flaring and cooling loops. The Giidel (1997) model is a modification 
of the solar nanoflare model proposed by Kopp & Poletto (1993). While the ASCA data cannot 
constrain many model parameters (number of loops, loop dimension, and mean magnetic field 
strength), the models can only reproduce the observed Em(T) provided that the power-law 
distribution of flare and microflare energies is a ~ 2. High S/N spectra of other active dwarfs are 
needed to test the viability of continuous flaring and to establish the importance of this mechanism 
as a function of rotation rate and effective temperature. 

In this paper, we present EUV and X-ray spectroscopy of HD 35850 = HR 1817 (F8 V). 
HD 35850 was detected serendipitously with EXOSAT (Cutispoto et al. 1991) and follow-up optical 
spectroscopy has established it as a nearby, single, solar-metallicity Pleiades-age (t ss 10 s yr), 
rapid rotator (Tagliaferri et al. 1994). We are particularly interested in HD 35850 because it is 
probably in the extreme state of magnetic activity for single, main-sequence F stars. There is 
considerable interest in probing the magnetic dynamo in F stars because they possess relatively 
shallow convection zones. Kim & Demarque estimate that, in a 10 7 yr-old 1.1-1.2M star like 
HD 35850, the convective turnover time is r c ~ 20 d. In a young 1.0M Q star like EK Dra, 
r c » 40 d. 
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2. Observations and Data Analysis 

HD 35850 was observed by the Extreme Ultraviolet Explorer from 1995 October 23 08:17:07 
UT to 1995 October 30 07:20:34 UT. The EUVE Deep Survey/Spectrometer consists of three 
aligned grazing-incidence telescopes. The telescope beams are intercepted by short-, medium-, and 
long-wavelength reflection gratings and detected with microchannel-plate detectors (SW, MW, 
and LW, respectively). The undispersed portions of the three beams are focussed onto the Deep 
Survey microchannel-plate detector, DS. This way, EUVE obtains simultaneous, time-resolved 
75-120 A broad-band photometry and 70-160 A, 170-370 A, and 300-525 A medium-resolution 
spectra (cf. Haisch, Bowyer, Sz Malina 1993). In support of this EUVE observing program, Call 
HK spectra were obtained on 1996 October 2-12. HD 35850 was observed with the ASCA SIS and 
GIS for ~ 18 ks on 1995 March 12. The reduction and analysis of these data are described below. 



2.1. EUVE Deep Survey Data 

The EUVE Deep Survey/Spectrometer records photon events whenever the source is visible 
by the satellite, including times when the satellite is passing through the South Atlantic Anomaly 
(SAA). These SAA passages are characterized by increased background levels which lead to 
significant and poorly characterized loss of telemetry, nicknamed "primbsching" . Correcting for 
earth occultations, the on-source time was ~ 203 ks. Correcting for instrument and telemetry 
dead times, the net exposure time was ~ 198 ks. 

The Deep Survey (DS) events were screened at various count-rate and primbsch-correction 
thresholds. Light curves were generated with bin times ranging from 100 to 2000 s. Source events 
were extracted from a 25 pixel radius circle centered on the PSF centroid. Background events were 
extracted from an annulus with inner and outer radii of 37 and 67 pixels, respectively. In Figure 1, 
we show the background-subtracted, dead-time and primbsch-corrected 75-120 A DS light curve 
for HD 35850 using 1000 s bins and including times when primbsching (monitor Det7Qldpc) was 
below 25%. 

The mean DS background-corrected count rate of HD 35850 is 0.20 counts s _1 . No large flares 
are evident in the DS light curve: the variability amplitude is approximately 100% and the most 
significant deviation from the mean count rate is a 5a peak near MJD 18.45. The light curve, 
however, shows small to moderate flares throughout the observation. At least 7 events have a 
peak count rate 3a above the apparent "quiescent" level of 0.20 counts s _1 . To test the hypothesis 
that the observed counts come from a non-variable source, we have determined the cumulative 
distribution function of the unbinned photon arrival times for the observed data and for simulated 
data from a constant, 0.2 counts s" 1 source. A two-sample Kolmogorov-Smirnov test statistic 
(Feigelson k, Babu 1992) was calculated using the observed and simulated distributions; the 
probability that the two data sets were drawn from the same parent population is low (P < 10 -9 ). 
We conclude that the DS light curve exhibits significant variability. 
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Fig. 1 — EUVE Deep Survey light curve (~ 950 s bins) of HD 35850 obtained 1995 October 23-30. 
The solid gray line indicates the mean DS count rate and the dashed gray line shows the estimated 
activity level during the ASCA observation in 1995 March. 
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Fig. 2. — Scargle discrete Fourier transform for the binned DS light curve in Fig. 1. The 
power spectrum indicates significant power at 0.0656 d, 0.98-d, and its aliases. These periods 
are associated with satellite SAA passages. The only peak which cannot be assigned to systematic 
effects occurs near 1.40 d. 
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We have used the binned DS light curve to compute the discrete Fourier transform (Scargle 
1989), shown in Figure 2. There is significant power around 0.98 d and its aliases. For HD 35850, 
vsini ~ 50 km s _1 and R+ ~ 1.18-Rq implies P rot /sini ~ 1.1 d, consistent with the 0.98-d 
period. However, Halpern & Marshall (1996) report that 0.98 d is a beat period associated with 
the satellite's passage through the SAA. For example, during the HD 35850 observation, the 
Det7Qldpc monitor shows strong periodic signals at 0.984 d and at 0.0656 d. Similar periods 
have been seen in other primbsch-corrected light curves (Marshall 1998, private communication). 
We thus conclude that some or all of the 0.98-d signal seen in the DS light curve results from 
SAA passages. We note that the only period seen in HD 35850's power spectrum that cannot be 
attributed to primbsching is at 1.40 d. Until further observations can be carried out, we tentatively 
identify 1.40 d as the rotation period of HD 35850 (see §2.2). 

2.2. McMath-Pierce Optical Spectroscopy 

Medium-resolution Ca II HK spectra were obtained with the Solar-Stellar Spectrograph at 
the National Solar Observatory's McMath-Pierce Telescope on 1996 October 2-12, ten days before 
the EUVE observation. The spectra around HK are typical of active, late-type dwarfs, with broad 
absorption troughs and bright, narrow emission cores. To calculate the equivalent width of the H 
and K emission, we measured the excess emission in the line cores, normalizing by the integrated 
continuum intensity in a 1 A box far from line center. 

In Figure 3 we plot the equivalent width of the Ca II H (dots) and K (asterisks) emission 
cores. While the equivalent widths are highly variable on short time scales, HD 35850 was only 
visible for a few hours each night. The limited sampling makes it difficult to estimate P rot with 
these data alone. In Figure 4, we have phase-folded the summed H and K equivalent widths using 
-frot = 1.40 d as derived from the DS photometry. 

2.3. EUVE SW and MW Spectrometer Data 

For the SW, MW, and LW spectroscopic data, the photon event data were screened to 
eliminate times of high background using the IRAF EUV package. Specifically, data were excluded 
when the total detector count rate exceeded 40 counts s -1 . The SW, MW. and LW photon event 
lists were used to create wavelength-calibrated FITS images and the SW, MW, and LW images 
were then analyzed within IDL. 

We note that during the HD 35850 observation, the SW microchannel-plate detector developed 
a hot spot 8-12 detector pixels from the dispersion axis, close to the location of the Fe XXII 114 A 
line. A light curve of the hot spot intensity was used to identify the times when the hot spot 
flickered; those times were removed from the photon event list. Some hot pixels were still visible 
in the SW image, so a 9 x 9 pixel box around the hot spot was set to the local background level. 
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Fig. 3. — McMath-Pierce Solar-Stellar Spectrograph Ca II HK light curves of HD 35850 obtained 
1995 October 2-12. For each exposure, we determined the excess emission in the line core in 1 A 
divided by the continuum flux in 1 A (i.e., the equivalent emission width). 
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Fig. 4. — Phase-folded light curve of the combined Ca II H and K equivalent widths for a period 
of 1.405 d and zero-point of JD = 2450000.445, the time of Ca II HK emission maximum. 
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Spectra have been extracted from the SW, MW, and LW images using a modified version of 
the IUE SIPS spectral extraction algorithm (Lenz & Ayres 1992). The optimal extraction routine 
is specifically tailored to low S/N spectra from photon-counting devices like the EUVE and IUE 
spectrometers. Briefly, a raw SW, M W, or LW 2048 x 2048 pixel image is trimmed and compressed 
to 1024 x 290 and 2-pixel Gaussian smoothed. Since the spectra are over-sampled (no less than 7 
detector pixels per resolution element), 2-pixel compression and smoothing improves S/N per bin 
without compromising spectral resolution. Using a predefined background region away from the 
source spectrum chosen to reduce curved edge effects caused by the EUVE reflection gratings, a 
background image is generated. A background-subtracted source image and corresponding error 
image are then created. The source image is used measure the cross-dispersion profile. To perform 
an optimal extraction, each row in the source image is weighted by the cross-dispersion profile and 
the source spectrum and corresponding error are extracted using a 7-pixel aperture. 

No emission lines longward of 365 A were detected in the LW count spectrum, presumably 
because of increased interstellar absorption at long wavelengths. Since the lines shortward of 
365 A are measured in the MW spectrum, we did not consider the LW spectrum in our analysis. 
In Figures 5 and 6 (bottom panels), we show the flux-calibrated SW and MW spectra. 

2.3.1. IDL Line and Continuum Fitting 

In order to measure line fluxes, we have constructed a line list using the emissivity lists of 
Brickhouse et al. (1995), Monsigniori Fossi & Landini (1993), and Mewe et al. (1985), in order 
of priority. Lines were grouped to match the SW and MW spectral resolution (0.5 and 1 A, 
respectively), allowing us to account for emissivity from blended lines. Based on the ionization 
states of the brightest identified lines, e.g., Fe XV, Fe XVI, and Fe XX, we made a starting guess 
at the plasma temperature (logT ss 6.8) and assumed a density n e = 10 11 cm~ 3 . The line 
emissivities, spectrometer effective area curves, and exposure time were used to provide a starting 
guess for each line strength. The observed spectrum was then fit with multiple Gaussians, varying 
wavelength (± one resolution element) and line strength. Lines with fewer than 2a counts were 
eliminated from the fit and the procedure was repeated. This way, line counts and upper limits 
were determined for all high-emissivity Fe IX to Fe XXIV lines. 

The SW spectrum shows a small but detectable thermal bremsstrahlung continuum. To 
determine the SW continuum luminosity from 85-135 A, the continuum bins were fit with two 
parameters (emission measure and temperature) using the continuum emissivity model of Mewe et 
al. (1985). The Hipparcos distance to HD 35850 is 26.8 pc (Perryman et al. 1997). For HD 35850, 
the 85-135 A continuum luminosity is L cont ~ 4.6 ± 0.4 x 10 28 ergs s -1 . 

In Figs. 5 and 6 (top and middle panels), we show the observed and fit count spectra with 
line identifications. Detected lines are indicated with dotted lines and upper limits are indicated 
with dashed lines. In Tables 1 and 2, we list the ion, possible blends, temperature of peak 
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Fig. 5. — EUVE SW spectrum of HD 35850. Top panel: observed count spectrum (0.135 A bins). 
Middle panel: detected (dotted lines) and undetected (dashed lines) SW lines plus continuum. 
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Fig. 6. — EUVE MW spectrum of HD 35850. Top panel: observed count spectrum (0.270 A bins). 
Middle panel: detected (dotted lines) and undetected (dashed lines) MW lines. Bottom panel: 
flux-calibrated MW spectrum. 
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emissivity, laboratory and measured wavelengths, a detection flag, measured line luminosities, 
and the signal-to-noise ratio for detected lines. The line luminosities in Tables 1 and 2 use 
A^h = 1-7 x 10 18 cm~ 2 (see §2.3.3). In all, 28 distinct lines are detected above 2a in the SW and 
MW spectra. 

2.3.2. Electron-Density and Column-Density Estimates 

The Fe XV and Fe XVI AA285, 335, 365 line ratios are relatively insensitive to density 
and provide the best estimate of interstellar Nr in the EUVE bandpass. We find 
N H = 1.9 ±0.4 x 10 18 cm" 2 . 

Two density-sensitive line ratios can be measured accurately from the SW spectrum of HD 
35850: Fe XXII AA114, 117 and Fe XXI AA102, 129. In Figure 7 we plot the observed SW spectrum 
(histogram) and fitted lines and continuum (solid line). The density-sensitive lines are indicated 
with dashed lines. The line Fe XXI line ratio is 0.22 ± 0.05. The predicted branching ratio (solid 
curve) and observed line ratio (cross with la error bar) are plotted versus density in Figure 8. 
Fig. 8 suggests that logn e < 11.6 cm -3 . Consequently, we have chosen logn e = 11.0 cm -3 for the 
IDL and SPEX emission-measure analyses. 

The Fe XXII line ratio is 0.47 =L 0.14, suggesting n c > 10 13 cm -3 . This ratio is unusually high 
and the implied density does not agree with the more reliable Fe XXI measurement. The detector 
hot spot (§ 2.3) may have contributed to the anomalous 114 A flux. We note that line ratio 
calculations may have systematic uncertainties of up to 50% (Brickhouse et al. 1995) and that the 
derived density limit is only approximate. 

2.3.3. IDL Emission- Measure Analysis 

The measured line luminosities and upper limits have been used to estimate the coronal 
emission-measure distribution, Em(T). Previous analyses of EUVE spectra (e.g., Giidel et al. 
1997; Mewe et al. 1996) have generally relied on global fits to estimate Em(T). For low S/N 
spectra like these, we prefer to use the bright lines in Tables 1 and 2 because (i) the emissivities 
of brighter lines have been fairly well established and (ii) the \ 2 statistic is not dominated by 
continuum bins. 

We have excluded the He II A304 line from our analysis. On the Sun, the He II Lya line 
is not a reliable temperature indicator: some, perhaps most, of the He II emission occurs in the 
chromosphere as a result of back- heating from the corona. We note also that around A255, the 
Fe XXIV line is blended with He II. The combined flux from both lines is listed as an upper limit 
in Table 2. For the remaining lines, our fitting procedure minimizes x 2 between the measured and 
predicted line luminosities by varying the amplitude and shape of Em(logT) (in AlogT = 0.1 




Fig. 7. — SW spectrum of HD 35850 in the 95-135 A region (histogram) and our fit to the lines 
and continuum (solid curve). The dashed lines indicate the wavelength of the density-sensitive 
Fe XXI and Fe XXII lines. 
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Fig. 8. — Fe XXI branching ratio as a function of logn c . The measured Fe XXI AA102, 129 
ratio and la error are plotted. Based on the predicted branching ratio, we derive an upper limit 
logn e < 11.6 cm -3 . 
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bins). 

We use exponential Chebyshev polynomials (Lemen et al. 1989) to describe Em(logT). 
In addition to detected emission lines, the undetected Fe X to Fe XIV lines have been used as 
additional constraints in the non-linear least-squares fitting. The SW continuum luminosity serves 
to constrain Fe abundance. 

Best-fit emission-measure distributions have been derived for a grid of column densities, Fe 
abundances, and Chebyshev-polynomial orders: 0.0 < Ah < 4.0 x 10 18 cm~ 2 , 0.1 < Z < 2.0, and 
2 < n < 11, where Z is the coronal Fe abundance relative to the solar photospheric value of Anders 
Sz Grevesse (1989). We obtain acceptable Xu values for all n > 6 and the resulting Em(T) are 
characterized by peaks at logT rs 6.8 and logT « 7.4. The best-fit column densities, abundances, 
and emission-measure distributions are not a strong function of n. As a result, we choose n = 6, 
the lowest polynomial order for which we obtain good fits. 

Fits using different order polynomials orders, column densities, and abundances yield 
qualitatively similar results: prominent peaks in the emission-measure distribution at logT ~ 6.8 
and logT ~ 7.4, although the sharpness of the high-T peak is not well constrained. The 
relative strength of the high-T component increases with increasing abundance. The best-fit 
emission- measure distributions all go to zero at very high temperatures (logT > 7.7). 

In Figure 9, we plot log A% 2 as a function of Ah and Z. The best fit (x 2 / ~ 1.07) is obtained for 
Nu = 1-7 x 10 18 cm* 2 and Z = 1.15. The 68% and 90% confidence contours are plotted for 9 free 
parameters (Z, Ah, and n = 0, . . . ,6). The EUVE line-to-continuum ratio yields Z = 1.15^']^. 
To obtain substantially sub-solar abundances (Z < 0.5), would require more than doubling the 
continuum level (dark solid line in Fig. 7). This is not compatible with the background-subtracted 
SW spectrum. We have experimented with various background-subtraction algorithms (smoothed, 
unsmoothed, polynomial fit): all produce bright lines and a weak continuum. 

We find adequate fits for a range of column densities and abundances, and, in particular, we 
find xl ~ 1-15 for Ah = 1.4 x 10 18 cm -2 and Z = 1. We adopt these values in Figs. 10-14. In 
Figure 10, the observed and predicted line luminosities for detected (filled circles) and undetected 
(open triangles) lines are compared, with a factor of two deviation indicated with dashed lines. 
Note the anomalously high Fe XXII line flux (see §2.3.2). In Figure 11, we plot Em(logT) per 
logT bin of 0.1 dex as derived from the IDL line analysis (dash-dotted line). 



2.3.4- SPEX Differential Emission-Measure Analysis 

In order to verify the results of the IDL line analysis, we performed a differential emission 
measure analysis of the optimally-extracted EUVE SW and MW spectra using the SPEX code 
(see Kaastra et al. 1992 for a detailed description of the DEM method and the SPEX Collisional 
Ionization Equilibrium model). For fixed values of Ah and Z, we derived Em(T) using the SPEX 
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Fig. 9. — Grid of log Ax 2 . The SW and MW lines and continuum luminosities were fit using 8th- 
order Chebyshev polynomials to describe Em(logT) at each point in the Z versus N# plane, where 
Z is the coronal Fe abundance relative to the solar photospheric values of Anders & Grevesse (1989). 
Dark pixels indicate regions of low \ 2 ■ The 68 and 90% confidence contours for 9 free parameters 
are shown. 
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Fig. 10. — Comparison of observed to predicted line luminosities for a range of Fe ionization states 
for Z = 1, iVn = 1.4 x 10 18 cm~ 2 , and logn e = 11.0 cm~ 3 . Downward-pointing triangles indicate 
upper limits for high-emissivity, undetected lines. Note the discrepant Fe XXII measurement. 
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Fig. 11. — Best-fit emission- measure distribution based on the IDL (dash-dotted line) and SPEX 
(solid line and error bars) analyses of the EUVE SW and MW lines and continuum. 
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regularization method. In Fig. 11, we plot Em(T) for N H = 1.4 X 10 18 cm" 2 and Z = 1 (solid 
line). The IDL and SPEX analyses yield similar emission-measure distributions, with 55-60% of 
the emission measure in the hotter component above log T = 7. 

2.4. ASCA SIS data 

An 18-ks ASCA exposure of HD 35850 obtained in 1995 has been analyzed by Tagliaferri et 
al. (1997). We use their reduced ASCA SISO spectrum to further constrain the EUVE results. 
Because the SIS cannot resolve individual emission lines from coronal sources, we cannot use 
the IDL line analysis method described in § 2.1.4. Consequently, we fitted the ASCA SISO 
spectrum using the SPEX DEM code. For the SISO fits, N u was fixed at 1.4 x 10 18 cm -2 . The 
best-fit abundance using the DEM collisional ionization equilibrium (CIE) model is Z ~ 0.5 with 
acceptable fits in the range 0.34-0.81. The ASCA upper bound on Z is plotted in Fig. 9 as a 
dashed line. We note that Z ~ 0.8 is marginally compatible with EUVE and ASCA. 

To illustrate, the ASCA SISO data are shown in Figure 12 (points) with the DEM model 
spectrum (histogram) for Z = 1. Discrepancies between the SIS data and the Z = 1 model are seen 
around 1.2 keV and 2.4 keV. These are the spectral features which normally drive the fit towards 
lower abundances. Brickhouse et al. (1997) have identified a complex of lines from 1.0-1.3 keV 
from highly excited (n > 5), highly ionized (Fe XVII to Fe XXV) states that are missing from the 
plasma codes in XSPEC and SPEX. The missing lines may partly explain the large discrepancy 
around 1.2 keV. 

In Figure 13, we show the corresponding ASCA SPEX Em(logT) for Z = 1. Like the EUVE 
emission-measure distributions, the ASCA Em(T) peaks at logT of 6.7 and 7.4. Finally, we note 
that the EUVE count rates are higher than expected from the SIS spectrum. Based on the SIS 
emission- measure analysis, the expected EUVE DS count rate is approximately 0.16 counts s _1 , 
indicated in Fig. 1 as a dashed line. HD 35850 appears to have been more active (by ~ 25%) 
during the EUVE observation. 

2.5. Comparison with Previous Results 

For HD 35850, we find that the EUVE line-to-continuum ratio indicates approximately solar 
coronal Fe abundance while a SPEX DEM analysis of the ASCA SISO spectrum is consistent 
with moderately sub-solar abundances. We note that HD 35850's photospheric Fe abundance is 
close to solar (Tagliaferri et al. 1994). Analyses of other active coronal sources often find coronal 
abundances far below the measured photospheric abundances (e.g., AB Dor, Mewe et al. 1996). 

The ASCA SIS and GIS spectra of HD 35850 were fit by Tagliaferri et al. (1997) using a 
number of multi-temperature plasma models in XSPEC (Arnaud 1996). For example, the best-fit 
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Fig. 12.— The 18 ks ASCA SISO spectrum of HD 35850 obtained 1995 March 12 (points and 
error bars). The best-fit SPEX DEM model (not shown) indicates sub-solar Fe abundances, 
0.34 <Z< 0.81. The model shown (solid histogram) assumes Z = 1, Nu = 1.4 x 10 18 cm -2 , 
and logn c = 11.0 cm -3 . The SISO spectrum is generally consistent with solar abundances aside 
from serious departures around 1.2 keV and 2.4 keV. 
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Fig. 13. — Emission-measure distribution based on the SPEX global-fitting of the ASCA SISO 
spectrum. The EUVE and ASCA data were not obtained simultaneously. Poor fits resulted when 
both data sets were modeled together. 
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SISO MEKAL parameters indicate a two-temperature corona {kT\ 0.6 keV and kT2 ~ 1.2 keV) 
with sub-solar abundances (0.12 < Z < 0.25). Tagliaferri et al. (1997) also performed a combined 
PSPC/SIS/GIS analysis using a 5-ks ROSAT PSPC observation of a field near HD 35850. The 
combined analysis involved comparing 2-T and 3-T MEKAL and RS models in XSPEC with solar 
and non-solar abundances and using different cross-detector normalization constraints (see their 
Table 5 and Figure 3). The XSPEC fits generally favor sub-solar abundance models; the best-fit 
MEKAL 3-T model yields Z = 0.34 ± 0.04, T t = 0.52 ± 0.07, T 2 = 0.78±g;J|, and T 3 = 1.9±^, 
with ~ 47% of the total emission measure in the coolest component. This cool component is the 
logT = 6.8 component seen in the EUVET5EM.. The two hottest components appear to represent 
the hotter DEM component. Our DEM analysis of the ASCA SIS spectra yields somewhat higher 
abundances (0.34 < Z < 0.81) than Tagliaferri et al. (1997) and somewhat lower Fe abundance 
than our EUVE analysis. 

Systematic differences between ASCA and EUVE analyses have been reported for other 
bright coronal sources. For example, Brickhouse et al. (1997) obtained simultaneous observations 
of Capella with ASCA and EUVE in 1996 March. As with previously obtained EUVE spectra of 
Capella, they find essentially solar photospheric Fe abundance. On the other hand, the plasma 
codes in XSPEC yield consistently poorer fits to the ASCA SIS spectra of Capella and require 
sub-solar abundances [Z ~ 0.7). 

Note: After the submission of this paper, a preprint of a paper to appear in Astronomy and 
Astrophysics by Mathioudakis & Mullan (1998) came to our attention which analyzes the EUVE 
observations of HR 1817=HD 35850. This paper and the Mathioudakis k, Mullan paper arrive 
at substantially different conclusions regarding the Fe abundance of HD 35850. Mathioudakis & 
Mullan (1998) visually compare the observed EUVE SW and MW spectra with synthetic spectra 
based on the 3-T MEKAL model suggested by Tagliaferri et al. (1997). They conclude that 
the EUVE spectra are consistent with Tagliaferri et al. (1997) and are inconsistent with any 
solar-abundance model, particularly because the Fe XIII— Fe XV lines in the MW spectrum are so 
weak. While it is difficult to visually compare spectra, their observed and synthetic MW spectra 
appear to be quite different: the predicted lines are weak and the the low-Z model appears to 
over-predict the MW continuum (see Figures 2 and 5 in Mathioudakis & Mullan 1998). Given 
the limitations of ASCA and EUVE, higher-resolution, higher-S/N AXAF HETG or XMM RGS 
spectra may be needed to determine HD 35850's coronal abundance. 

3. Microflaring 

The double-peaked EM distribution observed on three active, solar-type stars has been 
modeled by Giidel (1997) on the basis of a simplified stochastic flare model derived from the solar 
nanoflare model of Kopp & Poletto (1993). Rather than approaching the EM distribution problem 
hydrostatically, as is the case when using loop scaling laws, Kopp & Poletto (1993) use a simplified 
hydrodynamic model to treat a large number of flares. 
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The flare model reduces the loop hydrodynamics to a point model with a chromospheric 
energy sink. A series of energy pulses of finite duration is fed into the loop. This heating energy is 
lost by conduction into the chromosphere and by radiation into space. Without flaring, all loops 
are kept at the same equilibrium temperature dictated by hydrostatic loop scaling laws (Rosner, 
Tucker, & Vaiana 1978). Therefore, our model starts out at a lower threshold temperature 
and does not consider cooler loops. The salient feature of this model is its phenomenological 
similarity with more sophisticated hydrodynamic simulations in terms of the emission measure, 
temperature, and radiation history. The point model cannot treat the more complicated 1-D 
problem of radiation and conduction occurring between the corona, the transition region, and the 
chromosphere. The need to simulate a large number of flares requires such a simplified approach. 

Flares are ignited in each loop randomly distributed in time but satisfying a statistical number 
distribution in total flare energy E above a pre-defined threshold energy. I.e., dN/dE oc E~ a 
where N(E) is the number density of flares in the energy interval [E, E + dE]. On the basis of 
solar and stellar optical flare monitoring, 1.8 < a < 2.0 has been suggested by Hudson (1991). For 
consistency with Giidel (1997), we chose a threshold energy E > 10 27 ergs. In order to reduce the 
number of free parameters, all loops have the same height and vary only in their thickness between 
a minimum of 6.4 x 10 6 cm and a maximum of 4.3 x 10 9 cm to accommodate the range of flare 
energies. Each run produces a large number of small flares and only a small number of very large 
flares. Because no very large flares were observed by EUVE, the largest flares in the simulations 
were eliminated. 

The total EM distribution is determined by averaging (in time and space) the EM distribution 
of all loops. Similarly, the light curve is the superposition of coronal radiation leaving the star into 
2tt steradians, i.e., half the radiative losses are absorbed in the chromosphere, half are radiated 
into space. We have experimented with various parameters and find that our model reproduces 
the required coronal luminosity and EM distribution only if we assume compact loops, with a 
semi-length of 2.6 x 10 9 cm and essentially full surface coverage. The simulations provide fairly 
tight constraints on the power-lax index of flare energies (a ss 1.8). In this case, the equilibrium 
temperature is 6 x 10 6 K. 

Our best-fit microflare model for HD 35850 uses 1100 loops, distributed evenly over the star's 
surface, and accommodating 2.64 x 10 7 flares during a typical observation lasting 5 x 10 5 s. The 
resulting EM distribution is shown in Figure 14. The soft X-ray light curves from these simulations 
show the same 1-12 hour variability seen in ASCA SIS and EUVE DS light curves. Because the 
power-law distribution of flares favors small (microflare) events, the simulated light curves also 
show lower-amplitude, shorter-time scale variability seen in higher-S/N X-ray observations of dMe 
stars but not discernible for HD 35850. 

The model EM distribution exhibits the characteristic two-temperature structure with a 
broad intermediate minimum. However, in this model the EM structure is not just a result of 
radiative cooling as suggested by Gehrels & Williams (1993). It is a consequence of the balance 
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Fig. 14. — Emission-measure distribution produced by the microflaring model if we assume a = 1.8, 
compact loops, with a semi-length of 2.6 x 10 9 cm, and full surface coverage. In this case, the 
equilibrium temperature is T ss 6 x 10 6 K and a dominant peak occurs near T sa 26 x 10 6 K. 
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between flare heating, cooling by conduction, and cooling by radiation. Given our simplistic 
model, it is noteworthy that the two-temperature structure (or rather the broad minimum around 
10 MK) observed on many active normal stars can be explained, in part, by a hydrodynamic 
property of frequently-flaring loops. 

4. The Active Corona of HD 35850 

We summarize the combined EUVE and ASCA results for HD 35850 in Table 3. The 
ASCA SIS spectral analysis suggests sub-solar abundances and the EUVE line-to-continuum 
ratio indicates approximately solar photospheric abundances. Although Z f» 0.8 is marginally 
consistent with the EUVE and ASCA spectra, all the ASCA analyses yield systematically lower 
Fe abundance than the EUVE Fe line-to-continuum ratio. 

Assuming photospheric abundances, the corona is characterized by emission measure in a 
warm component from 5-8 MK and a hot component from 21-30 MK. Like other rapidly rotating, 
Pleiades-age dwarfs, HD 35850 appears to represent an activity extremum for main-sequence 
solar-type stars. HD 35850's inverse Rossby number is T c /P ro t ~ 20.7/1.40 ~ 15 and its X-ray 
surface flux is F x ~ 1.8 x 10 7 ergs s" 1 cm" 2 . For AB Dor (K0 V), F x ~ 1.7 x 10 7 ergs s" 1 cm" 2 
(Hempelmann et al. 1995) and r c /P rot > 120. For EK Dra (GO V), F x ~ 1.5 x 10 7 ergs s" 1 cm" 2 
and T c /P r ot ~ 15. AB Dor, HD 35850, and EK Dra appear to have saturated or nearly saturated 
X-ray activity and, despite having more rapid rotation and a deeper convection zone, AB Dor's 
X-ray surface flux is the same as HD 35850's. 

The EUVE DS light curve of HD 35850 shows about one moderate-amplitude (40% increase) 
flare per day. The X-ray and EUV variability and the presence of substantial emission measure 
above 20 MK suggests that some flare-like mechanism must be heating the corona. To test 
this hypothesis, we have modeled HD 35850's light curves and Em(T) using the hydrodynamic 
microflare point model of Giidel (1997). The simulations suggest a power-lax index of flare energies 
a w 1.8. 

On the Sun, transient brightenings seen in Yohkoh Soft X-ray Telescope images of solar active 
regions provide a direct measure of the microflaring amplitude and power-law index. Microflaring 
can provide, at most, 20% of the heating rate required to power the active-region corona (Shimizu 
1995). Ofman, Davila, & Shimizu (1996) have shown that the transient brightenings seen by 
Yohkoh may be a consequence of resonant absorption of global-mode Alfven waves in coronal 
loops, excited by random footpoint motions of these loops. A more recent analysis of Yohkoh 
data suggests that a steeper distribution of smaller flares, the so-called nanoflares, may also be 
present (Shimizu & Tsuneta 1997). Whether the microflares occur as a result of Alfven waves 
or magnetic reconnection, they cannot account for the radiative output of the Solar corona. For 
the rapidly rotating F and G dwarfs, however, the microflare models are able to reproduce the 
observed luminosity, EM distribution, and variability. 
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To better test microflaring on stars, accurate coronal density and abundance measurements 
as a function of temperature are needed to better constrain the emission-measure distribution. 
Photospheric magnetic field strength and filling fractions can further constrain loop models. Also, 
higher S/N, higher cadence hard and soft X-ray light curves will provide a statistically robust 
estimate of the distribution of flares. Planned observations of AB Dor with the Advanced X-ray 
Astrophysics Facility's High-Energy Transmission Grating Spectrometer will help resolve some of 
the density and abundance issues discussed in this paper and may lead to more quantitative tests 
of coronal heating models. 
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Table 1. HD 35850: 1995 October 23-30: EUVE SW Lines 
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Table 1 — Continued 
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Table 2. HD 35850: 1995 October 23-30: EUVE MW Lines 
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Table 3. HD 35850: Derived and Published Parameters 



Parameter Value Reference 
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